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ABSTRACT 

We determine the D/H ratio in the interstellar medium toward the DO white 
dwarf PG 0038+199 using spectra from the Far Ultraviolet Spectroscopic Explorer 
(FUSE), with ground-based support from Keck HIRES. We employ curve of 
growth, apparent optical depth and profile fitting techniques to measure column 
densities and limits of many other species (H2, Na I, C I, C II, C III, N I, N II, O I, 
Si II, P II, S III, Ar I and Fe II) which allow us to determine related ratios such 
as D/0, D/N and the H2 fraction. Our efforts are concentrated on measuring 
gas-phase D/H, which is key to understanding Galactic chemical evolution and 
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comparing it to predictions from Big Bang nucleosynthesis. We find column 
densities logAr(H I) = 20.41 ± 0.08, log A^(D I) = 15.75 ± 0.08 and logA^(H2) = 
19.33 ± 0.04, yielding a molecular hydrogen fraction of 0.14 ± 0.02 {2a errors), 
with an excitation temperature of 143 ±5 K. The high H I column density implies 
that PG 0038+199 lies outside of the Local Bubble; we estimate its distance to be 
297^^1;^ pc (la). [D I + HD]/[H I + 2H2] toward PG 0038+199 is 1.9li[J;i x 10"^ 
{2a). There is no evidence of component structure on the scale of Av > 8 km 
based on Na I, but there is marginal evidence for structure on smaller scales. The 
D/H value is high compared to the majority of recent D/H measurements, but 
consistent with the values for two other measurements at similar distances. D/0 
is in agreement with other distant measurements. The scatter in D/H values 
beyond ~ 100 pc remains a challenge for Galactic chemical evolution. 

Subject headings: ISM: evolution, stars: individual, white dwarfs, ultraviolet: 
ISM 



1. Introduction 

The ratio of D/H in the interstellar medium (ISM) offers a means to measure the effects 
of supernovae, stellar winds, infalling halo gas and ISM mixing on the chemical evolution 
of the Galaxy. The primordial D/H ratio is a firm prediction of the cosmic baryon density 
from Big Bang nucleosynthesis (Schramm & Turner 1998; Buries 2000; Buries et al. 2001). 
Measuring the D/H ratio in environments as close to primordial as possible in practice means 
making determinations in low metallicity, high redshift QSO absorbers. Kirkman et al. 
(2003), building on published work (O'Meara et al. 2001; Pettini & Bowen 2001; Levshakov 
et al. 2002) determined primordial D/H = (2.78^g 38) x 10~^ (lo")^ which is consistent with 
WMAP results and Big Bang nucleosynthesis (Romano et al. 2003; Cyburt et al. 2003). 
Closer to the Milky Way but presumably outside of the Galactic disk, the D/H ratio from 
the high velocity cloud Complex C is (2.2 ± 1.1) x 10"^ {2a) (Sembach et al. 2004), slightly 
lower than but still consistent with the primordial value. 

D/H ratios should be lower than the primordial value in regions where star formation 
has taken place, and thus, destruction of D I occurred in stellar interiors (astration). How- 
ever, fluctuations in the D/H ratio may occur due to a number of reasons. Typical processes 



^We will often distinguish between la and 2a errors in this work, given the variety of errors cited in the 
literature and some conventions observed in previous papers dealing with D/H. Unless otherwise stated, 
errors are la. 
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that could produce such D/H variations would include inhomogeneous ISM mixing by su- 
pernovae, stellar winds and the infall of relatively unastrated matter into the Galaxy. In 
addition, recent studies have been made of differential depiction of HD onto polycyclic aro- 
matic hydrocarbons (PAHs) or dust (Draine 2004; Lipshtat et al. 2004; Peeters et al. 2004), 
modulated over scales of hundreds of pc and several Myr (deAvillez & MacLow 2002). 

To probe the pattern of D I processing, measurements are being made of the gas-phase 
D/H ratio at a variety of distances. On the nearest scale, Linsky (1998) found D/H = 
(1.5 ±0.1) X 10^^ from 12 sight lines within the Local Interstellar Cloud (LlC). Further out 
in the Local Bubble, which is an irregular cavity filled with hot, low density, ionized gas 
with a radius ~ 65 — 250 pc (Sfeir et al. 1999), four sight lines are consistent with the LIC 
value (Oliveira et al. 2003; Moos et al. 2002, and references therein). Wood et al. (2004) 
determined D/H = (1.56 ± 0.04) x 10^^ {la standard deviation of the mean) for 22 sight 
hues in the Local Bubble compiled from the literature. Hebrard & Moos (2003) used 14 sight 
lines to find a value of D/H = (1.32 ± 0.08) x 10~^ from D/0 (building on the O abundance 
studies of Meyer, Jura, & Cardelli 1998), which is marginally in agreement with the LIC 
value.^ 

At distances between the Local Bubble and Complex C (0.1<d<3 kpc), however, D/H 
does show a dispersion of values. The lines of sight toward A Sco, S Ori, e Ori, LSS 1274, 
JL 9, HD 195965 and HD 191877 (York 1983; Jenkins et al. 1999; York & Rogerson 1976; 
Laurent et al. 1979; Wood et al. 2004; Hoopcs et al. 2003) were measured with Copernicus, 
the Interstellar Medium Absorption Profile Spectrograph (IMAPS) and FUSE to have low 
D/H values of 0.65 < (D/H)/10^^ < 1.00, with errors of <0.3 in D/H for each measurement. 
In sharp contrast, Feige 110 (Friedman et al. 2002) and 7^ Vel (Sonneborn et al. 2000) were 
measured with IMAPS and FUSE to have high values of D/H = (2.18l^J^) x 10"^ (90% 
confidence fimits) and D/H = (2.14 ±0.82) x 10~^ {2a), respectively. Although the number 
of D/H measurements outside the Local Bubble is not high enough to characterize the nature 
of the D/H dispersion, these two anomalously high values of D/H could indicate a significant 
inhomogeneity in the ISM. 

In principle, such inhomogeneities should be averaged out if one measures D /H over long 



common implicit assimiption is that the local disk value for D/H is approximately the same as the 

gas phase value. This misses any deuterium which is locked into grains or HD molecules, which may or 
may not be important. In one example, Jenkins et al. (1999) do not find HD to be a significant reservoir 
of D toward 5 Ori, which has a low D/H (0.74l^:}|, 90% confidence) and no HD lines in their IMAPS data 
(logAf(HD) < 12.8; Spitzer, Cochran, & Hirschfeld 1974). Lacour et al. (2004) also find that HD is not a 
significant reservoir of D within several kpc for a sample of seven stars. The issue is under current debate 
in the literature (Ferlet et al. 2000; Lipshtat et al. 2004). 
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enough sight hnes. However, confusion between velocity components at differing distances 
makes such measurements increasingly difficult for high H I column densities. Hebrard & 
Moos (2003) used D/0 and D/N measurements as proxies for D/H toward stars at distances 
d>100 pc, and concluded that D/H is generally lower outside of the Local Bubble. D/0 
and D /N measurements avoid the need to measure the H I column density, which can be 
difficult depending on whether the Lyman series falls largely on the fiat part of the curve 
of growth and/or whether data for Lya are available from HST or lUE. D I and O I also 
have column densities that only differ by ~ 2 orders of magnitude, which reduces systematic 
errors compared to D/H, where the difference in column densities is ~ 5 orders of magnitude. 
The D/0 ratio has its own complexities, though, due to a dearth of unsaturated O I lines 
for high column densities, and the possibility of uncertain O I oscillator strengths, whereas 
at the highest column densities, A^(H I) becomes easier to measure with the onset of Lyct 
and Ly/3 damping wings (Wood et al. 2004). 

In light of the D/H and D/0 measurements, it is therefore crucial to characterize the 
variations in D I abundance outside of the Local Bubble, both as a probe of astration, and 
as a means to constrain processes such as ISM mixing, infall of near primordial gas from the 
intergalactic medium and satellite galaxies and deuterium depletion. 

White dwarfs produce simple continua that act as ideal background sources for studying 
ISM absorption lines, and thus provide useful sight lines out to a few hundred pc. In 
this paper, we describe D/H measurements using FUSE observations of the hot DO white 
dwarf PG 0038+199. We present the target data in § 2. Details of observations and data 
reduction, including profile fitting (PF), curve-of-growth (COG) and apparent optical depth 
(AOD) techniques are in § 3. We describe our analysis in § 4, discuss the results in § 5 and 
summarize our conclusions in § 6. 

2. PG 0038+199 

The DO white dwarf PG 0038+199 was discovered via UV excess (Green et al. 1986). 
Wesemael et al. (1993) classified it as a hot DO star, and Dreizler & Werner (1996) performed 
a model atmosphere analysis. Dreizler et al. (1997) found log (H/He) ~ —0.7 from Keck 
HIRES data, though the H I emission line is very weak. Kohoutek (1997) listed the object 
as a possible post- planetary nebula. However, Werner et al. (1997) searched for an old 
planetary nebula in the region by direct Ha imaging, and found no evidence for one. A 
summary of stellar parameters is in Table 1. 

The distance to PG 0038+199 is estimated using the photometric parallax technique. 
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Using a magnitude V — 14.544 and color excess E(B-V)=0.037 (Table 1) and synthetic 
spectra computed in the optical, we calculate a distance of 2971Jq4 PC (with the errors 
dominated by the uncertainty in stellar gravity; see § 4.1 for model details). 

PG 0038+199 (Galactic £ = 119.79°, b = -42.66°) is outside of the Local Bubble, 
and is bright in the far UV based on lUE spectra. It is therefore a good candidate to 
help characterize the dispersion in D/H at large distances. Low dispersion spectra of 
PG 0038+199 exist, but are not suitable for measuring the H I column density from the Lya 
line. There are no other high resolution UV spectra of PG 0038+199 in the Multimission 
Archive at the Space Telescope Science Institute except for the observations presented 

here. 



3. Observations and Data Reduction 

3.1. Keck observations and reduction 

High resolution optical spectra were obtained with the HIRES echelle spectrograph 
(Vogt et al. 1994) on the Keck I telescope on 1996 July 24, using the red cross disperser to 
cover the wavelength region between 4260 A and 6700 A at a resolution 8.3 km s~^. The 
exposure time was 2000 sec. The standard data reduction as described by Zuckerman & 
Reid (1998) resulted in spectral orders that have a somewhat wavy continuum. To remove 
the echelle ripple, we used the spectrum of BD+33°2642 which was observed in the same 
night. Data were corrected to the vacuum heliocentric frame. Na I was used as an absolute 
velocity reference for the neutral gas in the FUSE data. 

3.2. FUSE observations 

PG 0038+199 was observed by FUSE for 13100 sec on 2000 July 27 (data set P1040201) 
through the large aperture (LWRS) in time-tag mode. For a description of the FUSE instru- 
ment see Moos et al. (2000) and Sahnow et al. (2000). Focal plane splits were not performed, 
because it was early in the FUSE mission and they were not routinely done at the time for 
exposures of that length. It was also observed for 6600 sec for program Ml 14, a calibration 
program for periodic monitoring of channel co-alignment, but those observations involved 
motions in the direction of dispersion, and therefore were not deemed scientifically useful for 
this work. 

Because the observations were done with the LWRS aperture, airglow is significant in 
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the troughs of H I Ly/3 and Ly7, with telluric O I A1027.5 A also affecting the Ly/3 damping 
wing. We excluded those regions from our analysis. 

3.3. FUSE data reduction 

The FUSE data were processed with the CalFUSE 2.2.2 pipeline. We summed the 
FUSE d&ta, interactively using software provided by S. Friedman, CORRCAL, registering each 
exposure based on cross-correlations around narrow absorption features near the middle of 
each channel segment. See Fig. 1 for the full spectrum. We then measured and corrected 
the a error array and wavelength calibration to prepare the data for profile fitting. See 
Appendix A for details. 

4. Analysis 

We prepared a grid of stellar models to provide a first estimate for the continuum for the 

FUSE data. We then used a suite of techniques, software and analyists to analyze the data, in 
an attempt to limit errors dependent on individual algorithms, programs or persons. In par- 
ticular, we employ curve-of-growth (COG), apparent optical depth (AOD) and profile fitting 
(PF) techniques at various points in the analysis, and present results both from the profile 
fitting codes vpfit (Webb 1987, see Bob Carswell's site http://www.ast.cam.ac.uk/~rfc/) 
and Owens (Lemoine et al. 2002) obtained by three different co-authors: C. Oliveira (CO) 
and G. Hebrard (GH) for Owens and G. WiUiger (GW) for vpfit. We stress that all 
analysis methods used the same data and error values. 

4.1. StellEir model grid 

To determine the effect of uncertainties in the continuum from the stellar model on 
column densities of ISM species, a grid of models was computed using TLUSTY (Hubeny & 
Lanz 1995) around the best-fit values for effective temperature and surface gravity T^ff — 
115000 K ± 11500 and log^ = 7.5 ± 0.3 (Dreizler & Werner 1996). We fixed [He/H] to 1000, 
because no photospheric H I absorption is seen in the FUSE data; the precise H abundance 
is difficult to constrain due to the strong ISM absorption and blending with He II absorption. 
The effect of any uncertainty on higher order D I lines for such a low H I abundance should 
be negligible. The models included nine energy levels for H I and twenty for He II, plus line 
blanketing caused by them. 
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Synthetic spectra containing only H and He for simplicity were computed over the FUSE 
range using the program SYNSPEC (Hubeny & Lanz 1995) and convolved with the instru- 
mental profile using the program R0TIN3 (provided with SYNSPEC). We find photospheric 
lines in the data from N IV A 923.7, 924.3, O VI A 1031.9, 1037.6, Si IV A 1128.3, Si V A 1117.8 
(blended with Fe VII), S VI A 1117.8 and Fe VII A 1117.6, 1163.9, 1164.2, 1166.2 (and pos- 
sibly 1146.5, 1155.0). We also computed a model with metals in addition to H and He to 
investigate the influence of metal lines on the Ly/3 profile (§4.6), with T — 115000 K and 
log^ = 7.5, [He/H]=5. (a conservative lower bound), [C/H]= 5 x 10-^, [N/H]= 5 x 10"=^ 
and [0/H]= 5 x 10~^, which is compatible with the Dreizler et al. (1997) upper limit on 
hydrogen abundance (not formally detected). 

For the purposes of fitting a stellar continuum to the data, the stellar radial velocity is 
estimated to be —52 km s~^, based on O VI photospheric absorption fines. We then fitted 
the model spectra described above to the FUSE data, using separate fiux normalizations 
for each of the eight FUSE detector/channel segments. The models are shown for the Ly/3 
region in Fig. 2. 

The stellar models were used as a first attempt at a continuum fit for absorption line 
analysis, except for LiFlB, which had to have a continuum fitted by hand due to the "worm" 
feature (Sahnow ct al. 2000), an optical anomaly due to the astigmatism in the system and 
a wire mesh in front of the detector. The stellar models were finally given a minor linear 
correction, determined by eye, to match the flux level in each channel segment. Further local 
zero-point corrections to the continuum will be allowed in the subsequent analysis. We see 
no indication of scattered light. 

4.2. Methods of line analysis 

Curve of growth: The COG technique is often used with FUSE data because it offers 
an independent technique to determine column densities and Doppler parameters for species 
with many different oscillator strengths, such as H2. See Oliveira et al. (2003) for details. 

Apparent optical depth: In the AOD technique, the column density is determined by 
directly integrating the apparent column density proflle over the velocity range of the ab- 
sorption proflle. See Savage & Sembach (1991) and Sembach & Savage (1992) for details. 
We use the AOD technique as a consistency check of the results obtained with PF and COG, 
and to determine lower limits on the column densities of species which only have saturated 
transitions in the FUSE bandpass. 

Profile fitting: If there are several components contributing to a given atomic transition. 
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profile fitting in principle yields much more information than simple curve of growth analy- 
ses. However, the FUSE resolution of 15 — 20 km s~^ is not high enough to resolve typical 
ISM lines, which generally arise from a number of narrow components over a small velocity 
interval. H I lines are particularly challenging, in that their intrinsic widths are high due to 
their low atomic mass. Nevertheless, for blended complexes of lines, one can obtain reason- 
ably accurate estimates for total column density (with a corresponding inflated composite 
Doppler parameter), provided that the distribution of intrinsic component optical depths 
and Doppler parameters is smooth (Jenkins 1986). We employ two profile fitting programs 
for this work: OwENS and VPFIT. Owens was developed by M. Lemoine and the French 
FUSE team at the Institut d'Astrophysique de Paris, and has been used for a number of 
D/H studies of Galactic targets involving FUSE data. See Lemoine et al. (2002); Hebrard et 
al. (2002); Hebrard & Moos (2003); Oliveira et al. (2003) for a fuller description. VPFIT has 
been used for many years, but its contribution to D/H work has so far concentrated on ex- 
tragalactic sight lines (e.g. Webb 1991; Carswell et al. 1996; Pettini & Bowen 2001; Crighton 
et al. 2003). We use both codes here as a check for the robustness of results between the 
two. 

VPFIT originally was developed for the analysis of ground-based echelle spectroscopy, but 
has also been used on various high resolution HST spectra. It generally works best with well- 
defined wavelength solutions, line spread functions and error arrays, which can be a challenge 
for FUSE data, and normally is not used to work with as many fitting windows as Owens. 
VPFIT also requires more steps than Owens to deconvolve the effects of turbulent and 
thermal velocity. Nevertheless, VPFiT is easily adaptable and well-supported. Its strength 
is that it can easily accommodate highly blended complexes, and can generate a grid 
relatively quickly compared to OWENS. For VPFIT results, our errors are calculated based 
on Ax^ = Xmin + 1 fo^^ and A^^ = Xmin + 4 for 2a. Like OwENS, VPFIT can allow 
a wavelength offset for each fitting region to accommodate uncertainties in the dispersion 
solution. Similarly, VPFiT can permit variations in continuum level. We leave such variations 
as zero point offsets in those windows which require them, because it is unlikely (and not 
desirable) that the continuum slope depart significantly from that of the stellar model over 
the short range (on the order of <1 A) of each fitting window. The rehance of VPFiT on 
a pre-determined continuum with zeroth order adjustments is one of its main differences 
from Owens, which uses polynomial approximations to model flips and dips in the local 
continuum and, if desired, to approximate absorption line wings. We employ VPFIT as a 
check for general OWENS results and the binning used in its analysis, and typically use 
VPFIT to solve for column densities and Doppler parameters for one or a small number of 
species at a time. For this reason, we generally leave the Doppler parameter free to vary 
by species, with the benefit that unphysical Doppler parameters may signal that the LSF or 
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blended component structure may be poorly understood in a fitting region. 

For the VPFIT analysis, the data were left unbinned, to permit the exclusion of the 
narrowest intervals which are suspected noise spikes and regions contaminated by airglow. 
For the AOD, COG and Owens analyses in this paper, the data were chosen to be binned 
by three pixels to increase the signal to noise ratio, as allowed by the oversampling of the 
LSF. For the VPFiT analysis, wc use negative a values as flags to omit data from profile 
fits, whereas for OwENS, bad data were flagged by high a values. Wavelength ranges for 
measurement intervals for the COG, AOD and the two PF methods were not required to be 
the same. 

All three methods described above provide formal statistical error estimates. Systematic 
errors, which were discussed in detail by Hcbrard et al. (2002), also have an important effect 
in overall error estimates, and will also be considered in this study. 

4.3. Na I Optical Data 

The Keck spectra of the Na I 5891,5897 lines are shown in Fig. 3. We used the PF and 
AOD techniques to determine the column density of Na I. Following the procedure outlined 
in Sembach k Savage (1992), we find equivalent widths Wa(A5891) = 173.72 ± 5.62 mA, 
and Wa(A5897) = 107.71 ± 4.00 mA. 

4.3.1. Profile Fitting 

Profile fitting of the Na I doublet with Owens was performed independently of the other 
atomic species, with one and two absorption components. The column densities derived 
with one and two components are consistent with each other: log A^(Na I) = 12.19 ± 0.04 
(2(t). The fits are shown in the bottom panels of Fig. 3. Using two components does not 
significantly improve the quality of the fit determined from the statistic, leading us to 
conclude that at the Keck spectrum's resolution of 8.3 km s~^ the Na I hues are consistent 
with a single absorption component centered at ~ —6 km s~^. (Radial velocities are 
given in the heliocentric frame). VPFiT gave a satisfactory fit for a single component at 
V = -5.1±0.3 km s-^ with Doppler parameter Z.^t^^ km s"^ and log A/"(Na I) = 12.25±0.04 
(2(7 errors). We adopt v = —5.1 km s~^ as an absolute velocity reference for the neutral gas 
in the FUSE data. We note that Na I is an imperfect tracer of H I (Cha et al. 2000), but 
have no other high resolution information for the ISM absorption velocity structure. 
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4-3.2. Apparent Optical Depth 

The top panel of Fig. 3 compares the apparent column density of the two components 
of the Na I doublet as a function of velocity. The two profiles are not consistent in the range 
— 12 < V < —1 km s~^, indicating that there is some unresolved saturated structure at these 
velocities. Disagreements between the two profiles for v < —20 km and v > +4 km 
are at the level of noise in these regions. We derive apparent column densities, log A^(Na l) = 
12.15 ± 0.02 for the weaker line (5897.56 A) and log A^(Na I) = 12.11 ± 0.02 for the stronger 
line (5891.58 A). Since the difference between the log of the two apparent column densities 
is less than 0.05 dex, we can use the procedure outlined by Savage & Sembach (1991) to 
correct for the small saturation effects in these transitions. A better estimate of the true 
column density is then achieved with logA^a = log A^(Na I)a5897 + ^ log A^(Na I)a5891-5897, 
where the first term in the equation is the column density derived from the weaker member 
of the doublet (Na I 5897 A) and the second term corresponds to the difference in column 
density derived from the two lines of the doublet. Using the equation above we obtain 
logA^(Nal) = 12.19 ± 0.04 (Icr), consistent with the profile fitting results and which we 
adopt. 

4.4. Moleculcir hydrogen 

The PG 0038+199 FUSE data have a rich H2 absorption spectrum, and many of the 
H2 features are saturated and/or blended. It is crucial, however, to have a good constraint 
on the various H2 transitions, because the lines are often blended with D I, H I and all other 
atomic ISM species. We used both curve-of-growth and profile fitting methods to constrain 
A^(H2). 

4.4-i- Curve of Growth 

For the COG method it was assumed that a single Doppler parameter, 6, describes all the 
rotational levels of H2 and HD. For both of these molecules, equivalent widths corresponding 
to a total of 91 different transitions (6, 7, 33, 25, and 14 for the J — 0, 1,2,3,4 rotational 
levels of H2 respectively; 6 for HD J = 0) were measured in all the FUSE channels where 
they are detected (Table 2). We did not use J = 5 for a COG analysis because it is too 
weak. The equivalent widths for each line, measured in different detector /channel segments, 
were then compared. All of them agree within the 2a errors, though we regard the results 
for J = and 1 with caution due to blending. These equivalent widths were then combined 
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through a weighted average, and a single-component COG was fitted to the 91 values. It is 
shown in Figure 4. The COG fit yields b = 3.11°;^ km s'^ 

4.4.2. Profile Fitting 

For profile fitting with Owens, CO used a single Gaussian line spread function (LSF) 
to describe the instrumental hne spread function, with a full width FWHM of 10.5 unbinned 
pixels. GH allowed the LSF to vary. Hebrard et al. (2002) found from OwENS profile fits 
for D/H work that Ax^ = 4 (or 2^) at the 2a level for errors in line fitting parameters such 
as column density, Ax^ = 9 (or 3^) for 3a, and similarly up to 7a. This Ax^ = k"^ for ka 
relation is used in our OwENS analysis to establish 2a errors based on averages from la to 
7a. The profile fitting of the H2 and HD lines was performed simultaneously with the other 
atomic species for which the column densities were being sought. We used 11, 18, 12, 8, 18, 
14 fines for the J = 0, 1, 2, 3, 4, 5 H2 rotational levels in all the FUSE chaimel segments where 
they were seen. H2 and HD were not restricted in radial velocity in comparison to atomic 
species nor, in the case of the analysis of GH, to each other, who found a radial velocity 
difference \v-h.2 — ^hdI = 4.6 km s"-*^ using this approach. We derive bB.2,HD ~ 2.7 km s~^. 

VPFIT is less well-suited than OwENS for fitting simultaneously the very large number 
of spectral regions containing H2 and HD lines in the PC 0038-1-199 data. We therefore used 
it to determine the wavelength-dependent LSF, based on H2 column densities and Doppler 
parameters from the COG and OwENS results (Appendix B). It is then possible to offer 
a consistency check and estimate of the profile fitting uncertainties with VPFiT, using a 
subsample of fitting regions from the LiFlA segment. We made trial fits fixing the Doppler 
parameter at 6 = 2.7, 3.1, 3.4 km s^^ based on the PF and COG estimates, and took the 
range of H2 column densities as a function of b for the various rotational levels as roughly 
2(7 systematic errors. The range in values from the COG and OwENS PF measurements 
were also treated as 2a contributions to the systematic error. We then compared the 2a 
statistical errors for the various PF and COG methods, and the variations in column density 
as a function of Doppler parameter and the variations in column density between the various 
PF and COG methods as approximate 2a errors in the approach taken by Wood et al. 
(2004), to form overall 2a estimates in the H2 and HD column densities. The exact H2 
column densities do not make a significant difference to atomic column densities except in 
the case of H I (§ 4.6). 

H2 and HD column densities are presented in Table 3. The quoted uncertainties arc 2a. 
Discrepancies in the COG and PF results make la uncertainties difficult to estimate, but 
we assume that la errors would be roughly half the 2a errors. The total H2 column density 
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along this line of sight is logN(H2) = 19.33 ± 0.04 {2a), and log(HD/H2) = -5.381^;^^ {2a). 
Column density results were insensitive to the exact LSF used. 

To estimate the H2 excitation temperature, we fitted the J — and 1 levels indepen- 
dently of the other levels, and find Tqi = 136 ± 21 K. The column densities for the various 
J levels divided by their relative statistical weights are shown in Fig. 5. We can also fit the 

H2 < J < 4 and HD J = column densities with a single temperature, T04 = 143 ± 5 K. 
We note that the < J < 4 levels are not usually consistent with a single temperature, and 
that T ~ 140 K is unusually high for Galactic H2 (Black & Dalgarno 1973; Snow et al. 2000; 
Rachford et al. 2001, 2002). Using the H I column density derived in § 4.6, we compute the 
fraction of H2 along the sight line, /(H2) = 2iV(H2)/(iV(H I) + 2iV(H2)) = 0.14 ± 0.02 {2a 
errors) . 

Finally, we checked the accuracy of our continuum placement by examining 33 HD 
lines in the six FUSE segments covering 915-1100 A. The HD lines are optically thin, and 
therefore quite sensitive to continuum level. VPFIT can calculate zero point continuum 
correction factors for each of the 33 fitting windows. A simultaneous fit to the HD hues 
yields an average correction factor of (cont) = 0.967 ± 0.044. The continuum we use is thus 
not systematically low or high, and adjustments have a scatter of ~4%. 

4.5. Column Densities of Metal Species 

The column densities presented in this work for the atomic species were determined 
with the PF, AOD and COG techniques. Wavelengths and fX values (and equivalent widths 
for Fe II) of the metal, H I and D I transitions used in this work are hsted in Table 4. In 
some cases, it was only possible to quote a lower limit via the AOD method for column 

densities due to saturated transitions in the FUSE bandpass. The COG was used only with 
Fe II, because this species is the only one which has enough non-saturated and non-blended 
transitions sufficient for this method. 

We included in our model fit for metal lines a molecular cloud with H2 and HD, as 
described in § 4.4, because many atomic transitions are blended with H2. The C I absorption 
will follow cool gas like that responsible for the H2 and Na I absorption rather than the warm 
gas that will be responsible for much of the H I, D I, N I and O I absorption. Therefore, for 
Owens PF results, for co-author CO, D I, N I, O I and Fe II were fitted in one component, 
and C I and H2 in another. Co-author GH did not include C I in his fits. Atomic column 
densities are listed in Table 5. 

An example fit to O I A974, which is the only available unsaturated O I line and the 
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entire basis of our O I column density measurement, is in Fig. 6. It is strongly blended with 
H2 J = 2 and 5. The resolution in SiC2A at 974 A is ~ 18 km s ^, which is among the 
best in our data. However, in SiClB it is ~ 30 km s~^, among the worst. The SiC2A data 
therefore provide the strongest constraints. Co-author GH only used the SiC2A data for his 
Owens fit, while co-authors CO and GW also included the SiClB data for their Owens and 
VPFIT (see below) calculations, respectively. All three co-authors found consistent results 
with each other. We regard A^(0 I) with caution, because of its dependence on one hne with 
good resolution only in the SiC2A data. We will discuss A^(0 I) further in § 5.1. 

For VPFIT PF results, wc used the wavelength-dependent LSF derived in Appendix B, 
and fitted each species independently. The Doppler parameters were left free to vary, as 
a check for LSF problems or evidence of substantial unresolved blending. A significant 
(> 2a) difference between VPFiT and Owens results only arose for N I, for which the best 
fit Doppler parameter would require h — Q.OIq^ km s~^, which is high compared to h from H2. 
If we fix the Doppler parameter to 6ni = 3.1 km s~^ (based on the value from H2), then we 
obtain A^(N I) from VPFIT, consistent with the fits from both co-authors who used OWENS. 
However, despite the simplicity of the Na I absorption, it is likely that the H I, O I and N I 
absorption (which will be from warm as well as cold gas) will have multiple components, 
in which case 6 = 9.0 km s~^ is not unreasonably high but is simply a consequence of the 
unresolved velocity structure. Likewise, assuming the h value from the H2 fit for N I may be 
completely inappropriate since to a large extent the H2 and N I absorption may be coming 
from entirely different places. 

Otherwise, the results are consistent between the one vpfit and two Owens calcu- 
lations, and in particular are independent of how C I is included in the fits. The column 
densities derived for the atomic species are presented in Table 5, and are a compromise of 
Owens and vpfit values. We also found a number of unidentified lines, fisted in Table 6, 
which may be either stellar or interstellar. 

4.6. H I and D I Column Densities 

The Lya transition provides the best determination of 7V(H I) for the highest column 
densities, based on damping wings. Unfortunately, only lUE low resolution data cover the 
Lya region for PG 0038-1-199, and the ISM Lya trough is filled in, presumably by geo-coronal 
emission. FUSE provides wavelength coverage from the Lyman limit up to Ly/5. We must 
therefore make an estimate of H I column density based on Ly/3, which is blended with a 
large number of H2 and O I transitions, but which offers promise due to strong damping 
wings. 



14 



4-6.1. HI and D I column densities 

H I: For the OwENS calculations, we fitted the Ljf3 profile separately from the other 
species. However, because H I and D I are blended with O I and H2, we used the column 
densities of the latter, determined in § 4.4.2-4.5, as fixed constraints. All four of the H I 
Ly/3 profiles obtained with the LiF and SiC channels were used simultaneously. The two LiF 
channels dominate the results because of their higher signal to noise ratio. 

The VPFIT analysis of the H I column density included a simultaneous fit to all segments 
covering Ly/5, Ly7, Ly8, Ly9 and LylO. However only Ly/? determines the column density, 
with the profile fit shown in Fig. 7. The higher order lines were included in the fit for 
consistency and to constrain the Doppler parameter, as will be discussed below. We did not 
include Ly7, 5, e and Ly-6 in the fits because they do not present strong enough damping 
wings to constrain the H I column density, and are too saturated to constrain the D I column 
density. As a check, we did obtain a statistically consistent fit for our adopted results (see 
below) for Ly7 and Ly(5 with VPFIT, after exercising care in excluding regions affected by 
geo-coronal emission, which we confirmed by separating our data into day and night subsets.^ 
The various H2 and HD column densities were fixed according to the values in Table 3, which 
resulted in a best fit of log iV(H I) = 20.41. Errors will be discussed in the next section. The 
effective H I Doppler parameter is 13.6 ± 0.2 km s~^ (2(7 errors), which is most likely an 
inflated value of h due to the presence of unresolved components (Jenkins 1986). 

D I: With Owens, the D I fines for Ly7, 8, 9, 10 were fitted, and the H I fine parameters 
were also permitted to vary in the fit. For GH, only the SiC2A data were used. A fit from 
Owens to Lyman 9 and 10 is shown in Fig. 8. The reduced value is 1.04, with 1442 degrees 
of freedom. CO included both SiClB and SiC2A in the calculation, and had = 1.65 for 
1878 degrees of freedom. Results between the two measurements were consistent, showing 
that including only the highest resolution data made no significant difference. 

With VPFIT, the reduced x^ value is 1.03 for 3635 degrees of freedom over Ly/3, 7,8, 9, 10, 
shown in Fig. 9. In case the saturated Ly/3 D I line skews the statistics, we also made a fit 
with only Ly7,8,9,10. Results were not changed significantly. 

Comparing the Owens and VPFiT results, we agree on a column density logA'"(D I) = 
15.75. Errors are discussed below. The Doppler parameter, assumed to be a composite from 
several unresolved components as for H I, is 6 = 11.9 ± 1.6 km s~^ (2(7 errors). 



■^Geocoronal emission for Ly7 has a weak red wing. Upon consultation with D. Sahnow and P. Feldman 
(2004, private communication), a similar effect is seen in LyJ and Lye emission. This is most likely an 
instrumental effect arising from geocoronal illumination of the entire LWRS aperture. 
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We checked the D I column density using the AOD and COG methods. For the AOD 

method we find logA^(Dl) = 15.68 ± 0.15 (la error in the mean). A COG study for 
the SiC2A data (which have higher resolution than SiClB) yields logAr(D I) = 15.801^;^^, 
b — lO.Tl^J km s~^. These values are consistent with our PF results. 



4-6.2. HI and D I Column density errors 

Hebrard et al. (2002) identified and discussed a number of systematic effects which can 
affect column density estimates. Following their example, we considered the potential sys- 
tematic effects from the following sources: 

Continuum uncertainty: Both VPFIT and Owens can explicitly include local normal- 
ization corrections in their calculations of statistical uncertainties. 

For the Owens analyses, the local continuum is free to vary for each spectral window, 
and is approximated by a polynomial of order ~ 2 — 4. Uncertainties from the continuum 
approximation are reflected in the statistics (e.g. Hebrard et al. 2002; Oliveira et al. 2003). 

For the VPFiT analysis, we allowed the continuum for each stellar model to have a zeroth 
order correction factor in each of the 12 fitting regions (4 for Ly/3 and 2 each for Ly7, 8, 9, 
10), with a mean value of 0.963 ± 0.041 (la). The correction factor is thus (just) consistent 
with unity, and is consistent with the continuum uncertainty from HD (§ 4.4.2). We note 
that the multiplicative corrections for the four Ljf3 regions are 0.974 ± 0.004, 0.930 ± 0.005, 
0.936±0.006 and 0.936±0.007 for the LiFlA, LiF2B, SiClA and SiC2B segments respectively, 
which may indicate a local systematic ovcrcstimation of the stellar model fiux. Continuum 
differences between the models throughout the data except in the Ly/3 region are <4 — 5% 
(§ 4.1), and are thus consistent with the local continuum variations in the Ly(3 region for our 
profile fits. Continuum differences in the Ly/3 region are discussed below. The stellar model 
including [He/H]=5 and metals produced smaller continuum normalization corrections at 
H I A923 and D I A922, because it included contributions from N IV. However, the model 
with metals did not produce a statistically acceptable overall fit in the Ljf3 region, in contrast 
to the metal-free models. In any case, the inclusion of metals in the stellar model had no 
significant effect on either H I or D I column densities. 

Stellar model: We compare results for the grid of stellar models described in § 4.1. The 
maximum deviation between the best fit and other models is ±10% . This occurs for the He- 
rich, metal-free models over 1025.2-1025.4 A (He II A1025), or for the model with [He/H]=5 
and metals at Ly/3. Typical deviations in the damping wings are on the order of < 5%. For 
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the Owens analysis, the Ly/3 region was fitted using the stellar model grid without allowing 
for continuum corrections, producing a systematic error of ~ ±0.05 dex {2a). For the VPFIT 
analysis, a similar calculation employing local continuum zeroth order corrections resulted 
in variations in column density of AlogA^(H I) = 0.02 and AlogA^(D I) = 0.002 (2(t). 
Line spread function: We varied the LSF by 2 km s^^ (~ Icr), which produced changes of 
AlogA^(Hl) = 0.005 and AlogiV(Dl) < 0.001. These differences are both significantly 
smaller than the statistical errors (see below) . We do not consider LSF uncertainties further 
in our source of error. 

H2 column densities: We varied the H2 column densities, trying the COG and PF results for 
log A^(H2) for the ensemble of J =0,1,2,3,4,5 and HD as rough 2a upper and lower log A^(H2) 
bounds, respectively. The resulting column density variations are AlogA'^(H I) = 0.03 and 
AlogA^(D I) = 0.01 (~ 2a). H2 column density uncertainty is a major systematic error for 
our H I and D I column densities. 

A grid in logA^(H I), logiV(D I) gives statistical errors of AlogA^(H I) = 0.02 and 
A log A'^(D I) = 0.04 (both 2a) . Combining the systematic and statistical errors in quadrature 
yields log Ar(H I) = 20.41 ± 0.08 and logAr(D 1) = 15.75 ± 0.08 {2a). 

We are encouraged that the results are so similar between the Owens and VPFiT anal- 
yses by three different people, and including four independent FUSE measurements. Fur- 
thermore, the results are robust, despite the use of binned vs. unbinncd data, the variation 
of continua, stellar models, LSF and H2 column densities. This line of sight is only the third 
D/H target for which only FUSE data are used to determine iV(H I), the first two being 
from Wood et al. (2004). 

4.6.3. Other estimates for HI column density 

The total H column density, IogiV(H I + 2H2), defined here as Ht, is 20.48 ± 0.07 {2a), 
which is in accord with that predicted using Na I (Ferlet, Vidal-Madjar & Cry 1985). They 
employed 78 stellar sight lines to determine logA^(Na I) = 1.04[logA^(H I + H2)] —9.09 (their 
equation 1), with a correlation coefficient of 0.85 and slope of 1.04±0.08. Using our measured 
column density of logiV(Na I) = 12.19 ± 0.04 {2a) yields logiV(HT) = 20.461};^^, with the 
uncertainty dictated by the slope error in the Ferlet, Vidal-Madjar & Cry relation. The color 
excess E{B -V) ^ 0.037 offers a prediction of logA^(Hr) = 20.3 with a scatter of ~ 50% 
(Savage & Mathis 1979). 

We also examined the H I 21 cm profile from the Leiden-Dwingeloo H I survey (Hart- 
mann & Burton 1997), and find logA'"(H I) = 20.36. This result is for a velocity resolution 
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1.03 km s ^, 0.07 K temperature sensitivity and 0.5° resolution on the sky centered at 
Galactic coordinates ^ = 120°, b = -42.5° (13 arcmin from PG 0038+199). See Fig. 10 
for a velocity profile of the H I 21 cm flux. We integrated at the position of PG 0038+199 
between —20 < f < 20 km s^^ to avoid a secondary peak of emission on the blue side. If 
we integrate over —100 < v < 100 km s~^, logA^(H I) = 20.47. The beam size is likely large 
compared to variations in H I, and there is no way to tell from our FUSE data how much H I 
lies behind PG 0038+199, which would increase the column density compared to what we 
could see in absorption with FUSE. Gas which is far behind PG 0038+199 would likely be 
very far from the disk, given its Galactic latitude. Despite these uncertainties, the closeness 
of the Leiden-Dwingeloo survey column density to that of our absorption measurement lends 
supporting evidence to our result, and implies that the H I is not very patchy in this part of 
the sky. In any case, there is no evidence for a column density as large as log A^(H I) = 20.57, 
which would be necessary to bring the D/H ratio (discussed in the next section) down to 
the Local Bubble value of 1.5 x 10~^. 



5. Results and Discussion 

Combining the results of our column densities for N I and O I with those for D I, HI, 
HD and H2, we obtain the relative column density ratios as listed in Table 7 and shown in 
Fig. 11. We consider the ratios between the column densities for D I, H I, O I, N I each in 
turn, beginning with ratios relative to hydrogen, and compare our results to those in the lit- 
erature. If we include H2 and HD contributions, we find log[iV(HT)] = 20.48 + 0.07 {2a), and 
log[iV(D I) +iV(HD)] = 15.76+0.08 {2a). We therefore quote various ratios including H2 and 
HD, as well as the traditional form D/H, plus the HD/H2 ratio. We approximate la errors 
for HD and H2 as half the 2a errors for contributions to [D I + HD]/[H I + 2H2], 0/[H+2H2] 
and N/[H+2H2]. We neglect HD for the total H contribution since (HD/[H I+2H2] < 10"^). 



5.1. Column density ratios 

[D I + HD]/[H I + 2H2]: We obtain a value of {l.dltoM) x 10"^ {2a errors), shown in 
Fig. 12, which is significantly high compared to D/H = (0.85 + 0.09) x 10~^ for sight hues 
with logA^(H I) > 20.5 and distances d > 500 pc (Wood et al. 2004). The possibility of 
a low D/H value toward distant sight lines was first proposed by Hebrard & Moos (2003), 
based on D/0 and D/N measurements, and previously published D/H values. Wood et al. 
corroborated the argument, compiling a large set of D/H measurements from the literature 
(though only four were for logA^(Hl) > 20.5). If we exclude the molecular component. 
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D/H = (2.19lo!5o) ^ (^c" errors) making an even greater discrepancy with similar H I 
column density sight lines. We believe that our [D I + HD]/[H I + 2H2] and D/H ratios for 
PG 0038+199 arc robust, based on measurements from three individuals using two different 
software packages and allowing for a variety of systematic errors. Our largest source of 
uncertainty is systematic errors in the H I column density, because we do not have access 
to Lya. We would require log A^(H I + H2) = 20.83 to have [D I + HD]/[H I + 2H2] = 0.85, 
which is many a from our result. Such a high value makes a poor fit to the data, and is 
shown in Fig. 13. 

[D I + HD]/0 I: Our result is (2.63^0:85) x lO'^ (2a). If the HD contribution is excluded, 
the ratio drops to: D/0 = (2.40+o;78) x lO'^ (2ct). For the sample logA'(D I) > 15.5 and 
distance d > 300 pc, our values are within the Hebrard & Moos range for (1.5<D/0<2.5) x 
10~^, shown in Fig. 14, though we note that there are only three values with that D 1 column 
density threshold in their sample. The A^(0 l) error dominates the uncertainty. 

0/[H+2H2]: We regard our O I column density measurement with extreme caution, be- 
cause it depends entirely on one line at 974 A strongly blended with two H2 lines, dominated 
by the SiC2A data. We compare our results with the 0/H determination of Meyer, Jura, 
& Cardelli (1998), whose sample resembles our sight line in H I column density.^ Our value 

of O/P+2H2] = {7.76tlil) X 10"^ (2a errors) is > 2a higher than the local ISM value of 
0/H = (3.43 ±0.15) X 10^"^ (la) from Meyer et al. using the updated O I oscillator strength 
of Welty et al. (1999). If a H I column density error were the sole cause of the 0/H discrep- 
ancy, we would require logA^(H I) = 20.82 to bring 0/H to 3.58 x 10~^, or within the la 
error of the revised Meyer, Jura, & Cardelli value. However, such a high value for iV(H I) is 
highly unlikely, given our determinations from the FUSE data and in light of other A^(H I) 
estimation methods discussed in § 4.6.3. 

Another source of 0/H uncertainty depends on the accuracy of the O I 974 A line 
oscillator strength, which is theoretically calculated and is estimated to be good to 25% 
(Biemont & Zeippen 1992; Wiese, Fuhr & Deters 1996). In addition to PG 0038+199, 
other objects in the literature for which A'"(0 I) depends only on O I 974 are JL9, LSS 1274 
(Hebrard & Moos 2003; Wood et al. 2004) and HD 90087 (Hebrard et al. 2004). ^ They, 
too, exhibit high values of O/H, though the local ISM value hes close to or within their 2a 



''The samples of Andre et al. (2003) and Cartledge et al. (2004) have significantly higher characteristic 
H I and H2 column densities than PG 0038+199. 

'"'HD 195965 (Hoopes et al. 2003) also has high 0/H, but the H I and O I column densities have two 
independent measurements each from HST, FUSE and lUE, and the star's position in the Cygnus 0B7 
association allows for the possibility of containment within a metal-enriched cloud. 
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formal errors. If the O I 974 oscillator strength were to increase by 15% (0.06 dex), then the 
nominal O I column density toward PG 0038+199 would drop to logA^(0 I) = 17.31, with 
O/P+2H2] = 6.8 X 10"^, putting the local ISM value within the 2a error for PG 0038+199. 
The nominal D/0 value would increase from to 2.40 x 10~^ to 2.57 x 10^^, which still 
lies within the scatter of high column density data points from Hebrard & Moos (2003) 
(see above). The O I 974 oscillator strength appears consistent with that of O I 1356 in 
the case of HD 195965 (Hoopes et al. 2003; Hebrard et al. 2004), with the Hoopes et al. 
(2003) uncertainty in logA^(Ol) being +0.04,-0.06 dex - close to the value we require 
if the discrepant 0/H value toward PG 0038+199 arises because of an oscillator strength 
difference. Our D/H and D/0 results are consistent with a larger variation in D/H than 
in D/0 (Hebrard & Moos 2003). However, we find 0/H toward PG 0038+199 anomalous, 
because the variations in 0/H found by Meyer et al. were small, and we cannot identify a 
specific reason why the measurement would be in error. If we exclude the H2 contribution, the 
0/H discrepancy is even larger: 0/H = (8.5ll38g) x 10~"^ (So"). Observations of additional 
lines of sight will show whether the combination of D I, H I and O I column densities found 
toward PG 0038+199 produce an outher, are consistent with other sight lines or are subject 
to an unidentified error. 

D/N: We find [D I + HD]/N I = (2.341^;^^) x 10"^ {2a), which is high compared to the 
Hebrard & Moos (2003) values of (1.0<D/N<1.5) x lO'S shown in Fig. 15, for logiV(D I) > 
15.5. Omitting the HD fraction results in D/N = (2.29lo:??) x 10~^ (2a). 

N/P+2H2]: We find values of (0.81^°;^^) x 10"^ and N/H = (0.951°;^^) x 10"^ (all 2a), 
which are consistent within la errors of the ISM value of N/H = (0.75 + 0.04) x 10~^ (Meyer, 
CardeUi, & Sofia 1997). 

N/0: Our value of N/0 = (1.12^^;^^) x 10"^ {2a) is low compared to the ratio (2.35 ± 
0.15) X 10-^ for N/H and 0/H in the ISM from Meyer, CardeUi, & Sofia (1997) and Meyer, 
Jura, & CardeUi (1998), due to the high 0/H ratio we find. See the discussion above about 
the reliability of the O I measurement. 

HD/2H2: The deuterium fraction in molecular form is (2.08^025) x 10^^ {2a errors), 
and the column density ratio is HD/H2 = (4.17^^^;^^) x 10"^ {2a). The HD and H2 column 
densities place it within the range of measurements at similar A^(H2) from Spitzer, Cochran, 
& Hirschfeld (1974) and Savage et al. (1977), a factor of ~ 3 above the expected ratio from 
charge exchange reactions between H and D (Liszt 2003), but consistent with the scatter in 
values. 
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5.2. Compcirison with other high D/H Unes of sight 

There are only two other noteworthy detections of D/H > 2.1 x 10^^: 7^ Vel (Son- 
neborn et aL 2000) and Feige 110 (Friedman et al. 2002).^ 7^ Vel exhibits both high and 
low ionization absorbers in a 7-member complex (3 H II, 4 H I). Its H I absorption spans 
a velocity interval of Av — 32 km s~^, with a very low H2 fraction of 1.4 x 10~^, and 
logA^(H I) = 19.71 ± 0.03 (1(7 errors), A^(D I) = (l.Utoll) x 10^^ (90% confidence limits) 
and D/N = (2.7 ± 0.4) x 10"^ at 90% confidence or 1.65(7 significance (Cha et al. 2000; 
Fitzpatrick & Spitzer 1994). Although the D/N ratios are consistent within the errors for 
7^ Vel and PG 0038+199, in contrast, PG 0038+199 does not exhibit numerous high ioniza- 
tion species (though it does show S III), and has a much higher H2 fraction and H I column 
density. We cannot compare D/0 because the O I column density is difficult to measure for 
both 7^ Vel and PG 0038+199. Feige 110 has only a slightly smaller H I column density than 
PG 0038+199 (logAr(H I) = 20.14l°i^^ 2a errors), and D/H = (2.14 + 0.82) x 10"^ {2a). It 
is possible that PG 0038+199 is slightly higher out of the Galactic plane, and is samphng 
richer, infalling gas, one of several explanations which we consider in the next section. 

Another way to characterize the nature of the PG 0038+199 sight line versus the 7^ Vel 
sight line is in terms of the iron abundance. Jenkins, Savage, & Spitzer (1986) have examined 
the iron abundances determined from N(Fe II) measured by Copernicus as a function of the 
average density {Nh). (Nh) is defined as iV(H) = A^(H I) + 2A^(H2) divided by the distance 
to the star. For PG 0038+199, (Nh) = 0.3+o:i cm'^ (la) and the relative iron abundance 
assumed to be log(iV(Fe II)/A^(H)) = —6.06. This iron abundance is higher than average, but 
consistent within the uncertainties with the mid-density range results presented by Jenkins 
and coworkers. On the other hand, 7^ Vel has an iron abundance more than twice as great 
as PG 0038+199 (Jenkins, Savage, & Spitzer 1986; Fitzpatrick & Spitzer 1994) and {Nh) 
is much lower, 0.06 cm~^. In the model discussed by Spitzer (1985) and used by Jenkins, 
Savage, & Spitzer to analyze the Fe abundance, the two sight lines are quite different. The 
7^ Vel sight line samples only relatively warm diffuse gas with a density of ~ 0.1 cm~^, while 
the PG 0038+199 sight line samples both warm gas and cold higher density fluctuations with 
a density of ~ 0.7 cm^'^ (Spitzer). Jenkins, Savage, & Spitzer assumed that the properties 
of the grains vary between the two densities and showed by fitting their data set that the 
relative iron abundance would decrease a factor of four between the extreme cases of only hot 
gas and only cold gas. However, in this case, although the iron abundance changes by more 
than a factor of two, as expected, these two sight lines do not show an analagous change in 



Cru has high D/H but very large errors: 2.5tQ;g x 10 ^ (presumed la errors). It has no reported O I 
or N I measurements, and has by far the lowest H I column density among the high D/H sight lines, with 
logiV(H I) = 19.6 ± 0.1 (presumed la, York & Rogerson 1976). 
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the abundance of deuterium. The values of (Nh) and the Fe abundances are quite different 
for PG 0038+199 versus 7^ Vel, indicating that the grain properties along the two sight lines 
are likely to be different, but there is no evidence for a similar variation in the deuterium 
abundance. 



5.3. Variations in D/H: possible causes 

A number of different potential causes for variations in D/H have been discussed in 
the literature. All causes that have been proposed could produce inhomogeneities in D/H, 
provided that the ISM mixing time is sufficiently large. Wood et al. (2004) suggested that 
D/H variations occur on a scale where ISM mixing produces inhomogeneities on scales larger 
than the Local Bubble, which has a well-established D/H ratio, yet small enough so that the 
sight lines to the most distant targets investigated for D/H pass through enough inhomoge- 
neous zones for D/H variations to average out to their observed low value of 8.5 x 10~^. 
The strongest challenge to the Wood et al. model would be a sight line with A^(H I) > 20.5 
(or in the Local Bubble) with a high D/H value. PG 0038+199 is the highest H I column 
density sight line with a high D/H value, just below the logA'"(H I) = 20.5 upper hmit for 
the meso-mixing scale of Wood et al. We now briefly discuss some of the proposed causes of 
D/H variability below, but see Lemoine et al. (1999) and Draine (2004) for details. 

a. Inhomogeneous Galactic infall: Material which has undergone little astration could 
fall inhomogeneously onto the Galactic disk and produce the observed variation in D/H 
distribution. One example is from Chiappini, Rcnda, & Matteucci (2002), which predicts a 
gradient in D/H as a function of Galactic radius. More data are needed to determine the 
validity of the model (e.g. Hebrard & Moos 2003). 

b. Astration variations: Variations in the amount of stellar processing experienced by 
different regions of the ISM could result in D/H variability. An anticorrelation between D/H 
and 0/H has been suggested, on the basis that gas is enriched in O as D is destroyed in 
stars (e.g. Steigman 2003). However, the suggestion was tentative because the D/H sample 
size at the time was noted to be limited. PG 0038+199, with its high values of both D/H 
and (if confirmed) 0/H, clearly would be an exception to such an anticorrelation. We share 
Steigman 's view that more sight lines need to be analyzed to explore the topic further. 

c. Variable deuterium depletion: Variations of D/H could be caused by preferential 
depletion of D onto dust grains relative to H. Draine (2004) noted that such a process could 
account for the variation in D/H in the ISM while allowing for the relative homogeneity of 
abundances for N and O. Lipshtat et al. (2004) showed that it is possible for dust to enhance 
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the production rates of HD and D2 vs. H2. Observational evidence comes from Peeters et al. 
(2004), who found infrared spectral evidence of deuterated polycyclic aromatic hydrocarbons 
(PAHs), and pointed out the large deuterium excess in primitive carbonaceous meteorites 
and interplanetary dust particles (e.g. Messenger 2000). Draine also suggested that high 
velocity shocks could also work in the opposite sense to destroy grains and increase D/H. 

Dust depletion processes occur in cold regions of the ISM where H2 is the dominant 
reservoir of H. Preferential depletion of D onto dust requires particularly low temperatures 
(Draine 2004), so one might expect lines of sight with low H2 temperatures to have low D/H 
ratios. Typical ISM H2 temperatures have been measured to have mean values of Tqi ~ 70 K 
and rms scatter of ~ 15 K for over 80 hues of sight (Savage et al. 1977; Rachford et al. 2002). 
Wood et al. (2004) noted similar temperatures toward the high H I column density sight hues 
JL 9 and LSS 1274, and proposed that hues of sight which were recently shocked should have 
relatively high D/H. H2 temperatures could be increased in star formation regions, such as 
near a Cru and 7^ Vel, the latter being near the Vela supernova remnant. As gas cools, 
deuterium could be preferentially depleted onto PAHs or grains (Draine 2004; Peeters et al. 
2004; Lipshtat et al. 2004). The H2 temperature of Tqi = 136 ± 21 K toward PG 0038+199 
is unusually high compared to the mean ISM value, so this is consistent with a high D/H 
value in the depletion model. Quantitatively, such depletion could reduce D/H by 1 x 10^^, 
which is consistent within the errors with the difference in D/H between the PG 0038+199 
sight hue and sight lines with logiV(H I) > 20.5. 

However, we find no evidence of supernova remnants near PG 0038+199 (Stephenson & 
Green 2002), nor any sign of nearby anomalous X-ray emission (Snowden et al. 1995) or Ha 
emission (Haffner et al. 2003). The H2 temperature of PG 0038+199 could be raised by the 
faint remnants of a planetary nebula, though. Two other similarly hot DO stars have Ha 
detections of huge, old, faint PN: PG 0109+111 (Werner et al. 1997, a search which found 
no old PN around PG 0038+199) and PG 1034+001 (Rauch, Kerber, & Pauh 2004). 

If a correlation exists between H2 temperature and deuterium depletion, it would have 
to explain sight lines with low D/H and high H2 rotation temperatures. 5 Ori and e Ori 
(Jenkins et al. 2000) exhibit such high rotation temperatures with absorption components 
having T03 ~ 190 — 850 K. (The temperatures for Tqi are similar to T03 for those sight 
lines. Reasons for the similarity are discussed in Jenkins et al., though we recognize that 
Tqi is usually assumed to be the best estimate of the actual thermal temperature.) The H I 
column densities toward 6 Ori and e Ori are relatively high (logA^(H I) = 20.193 ± 0.025, 
log7V(H I) = 20.40 + 0.08, la errors), yet the D/H ratios are low {OMtoH, 0.65 + 0.3 x 10"^ 
respectively, la errors). Their H2 fractions are low, however (7 x 10~^, 1.3 x 10~^), so a 
different physical mechanism from that toward PG 0038+199 may be at work. In addition. 
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Spitzer, Cochran, & Hirschfeld (1974) found the rotation temperatures of 7^ Vel and ( Pup 
to be To3 = 964 ± 119 K and T04 = 1140 ± 154 K, (both la) respectively. (Spitzer, Cochran, 
& Hirschfeld state that their excitation temperatures for T03 and T04 fit all the observed 
values of A^(H2) ^or the J levels included in the calculation; we estimate Tqi from their H2 
column densities for J = and 1 and their equation 1 to be ~ 300 ±150 K and ~ 350 ±100 K 
respectively, which are still high.) D/H toward 7^ Vel is high, as discussed in the previous 
section, but toward C Pup is {lA2toll) x 10"^ (90% confidence, Sonneborn et al. 2000), 
which, given our uncertainties, is marginally below the value for PC 0038+199. A larger 
data set is necessary to investigate the relationship of H2 temperature and D/H ratio in 
further detail. Additionally, if this depletion scenario is correct, then D/H may be correlated 
with the abundance of refractory elements such as Fe. However, as discussed in § 5.2, there 
is no convincing evidence of such a relationship for PC 0038+199. 

It is clear that variation in D/H is significant outside of the Local Bubble, and its cause 
remains an open question. Our understanding of variations in the D/H will improve as 
theory evolves and further measurements of D/H, H2 and other ISM species are made at 
large distances in the Galactic disk and halo. 

6. Conclusions 

We have observed the He-rich white dwarf PC 0038+199 with FUSE and Keck+HIRES 
and found the following. 

1. Keck HIRES data hint at unresolved structure on the order of < 8 km based on 

Na I. 

2. We have compared column densities for a number of ISM species based on a variety of 
measurements using profile fits from OWENS and VPFIT, curve-of-growth and apparent opti- 
cal depth determinations, and found broadly consistent results between the various methods. 

3. The H I, D I and H2 column densities are log A^(H I) = 20.41 ± 0.08, log A^(D I) = 
15.75±0.08, log A^(H2) = 19.33±0.04 with a H2 fraction /jj^ = 0.14±0.02 {2a errors quoted 
throughout the conclusions, unless otherwise noted). The measurement of A^(H I) for a D/H 
target is only the third based solely on FUSE data. The H2 fraction is high enough that we 
include effects for it and HD in our column density ratios. 

4. [D I + HD]/[H I + 2H2] = {IMtolD X 10"^ is unusually high for a Galactic value. 
Variations in D/H outside the Local Bubble are significant. N I/(H I + 2H2) = (0.8llo:29) x 
10-^ is consistent with the mean ISM value (Meyer, CardeUi, & Sofia 1997). [D I + HD]/N I = 
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(2.341^:1) X 10-^ is high for log Ar(D I) > 15.5 (Hebrard & Moos 2003, though their sample 
size is small). 0/[H+2H2] = (7.76+|^^) x 10"^ is unusually high, with logA^(0 I) depending 
on a single line measurement. However, D/0 = (2. 40^*^0 78) ^ 10~^ is consistent with three 
other distant sight lines (Hebrard & Moos 2003). Further observations should determine 
whether the high 0/H value remains an anomaly, arises as a result of some uncertainty (e.g. 
O I 974 oscillator strength) or can in fact be statistically expected. 

5. The HD/H2 ratio is (4.17^049) ^ 10^®, consistent with the scatter of previously 
measured values for logA^(H2) ~ 19.3 (Liszt 2003). 

6. A number of mechanisms which produce variation in D/H remain plausible. The H2 

excitation temperature of 143 ± 5 K (la) may indicate past star formation in the vicinity of 
PG 0038+199, liberating D which may otherwise be depleted onto grains, thereby potentially 
explaining the high D/H value. However, there are low D/H sight lines with even higher 
excitation temperatures. More sight lines outside the Local Bubble need to be observed to 
characterize the D/H distribution and physical processes associated with it. 
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work has been done using the profile fitting procedures OWENS, developed by M. Lemoine 
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native Hawaiians for enabling observations from Mauna Kea. 



A. Details on the reduction of FUSE data 

We present details on the reduction of FUSE data, as the use of vpfit creates a need 
for accurate a errors, a desire for the most accurate wavelength calibration possible and 
requirement for knowledge of the LSF as a function of segment and wavelength. 
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A.l. Correction of a values 

We established a minimum value for the a array based on the distribution of a values 
for each segment. Individual pixels with values of a which are unrealistically low make it 
difficult or impossible to obtain acceptable statistical results for profile fitting. The mean 
flux and signal to noise ratio per pixel for the eight segments before modification of the a 
values, for A > 915 A and where cr > 0, ranges over 1.4 — 1.8 x 10~^^ erg cm~^ s~^ A"^ 
and 5.0 — 18.2.^ We plotted histograms of the unmodifled a values, and noted narrow, 
secondary peaks of low a values significantly (4-5 full width half maxima) below the main 
peak of the a distribution. There were a number of instances with (7 = 0. We therefore 
set a minimum a value for each FUSE detector segment. The fioors in the a distribution 
range from 1.4 — 5.5 x 10^^^ erg cm~^ s~^ A~^, which affects pixel fractions of 0.4-1.5% per 
detector segment. An example is shown in Fig. 16. 

We then calculated an empirical correction factor to the la error arrays propagated by 
the CalFUSE pipeline and affected by any rebinning of the data, fixed pattern noise etc. 
The correction is necessary to ensure proper probabilities for the profile fitting, and is 
recommended in VPFIT documentation. We first divided the data in each detector segment 
by the continuum to normalize the spectrum. We then binned the data for each segment, 
starting with increments of 5 pixels, divided the root mean square (rms) by the mean of the 
la error array (a) for all bins in a segment, and examined the distribution of rms /a ratios 
as a function of bin size to find a characteristic ratio. The process was repeated for bins of 
6,7,8, ...,100 pixels, in increments of one pixel. 

The mode of the distribution is the most stable characteristic quantity of the distribution 
of rms /a ratios, because the mean and median can be skewed by bins containing large 
absorption features. We therefore plotted the mode of rms/ (a) as a function of bin size to 
look for a trend in the values. A resolution element is ~ 10 unbinned pixels, therefore the 
mode grows quickly with bin size up to ~ 10 — 15 pixels for small bins, then stabilizes as 
the bin size increases to 20-40 pixels, which is characteristic of the size of spectral intervals 
between absorption lines. (There are deemed to be no significant intrinsic emission lines 
in the stellar spectrum.) At the bin size >40 pixels, the ratio increases again because 
the wavelength range covered becomes characteristically larger than the interval between 
absorption features. There was always a plateau in the value of the mode at ~ 20 — 40 pixels 
in bin width. We took the value of the plateau to be the a correction factor. The correction 
factors (plateaus) are 1.4, 1.1, 1.3, 1.3, 1.1, 1.1, 1.2, 1.1 for LiFlA, LiFlB, LiF2A, LiF2B, 
SiClA, SiClB, SiC2A, SiC2B respectively. 



''Unbinned pixels, 0.0066 A/pixel, are used throughout this analysis unless otherwise noted. 
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The large correction factor for LiFlA is likely due to fixed pattern noise. Because focal 
plane splits (which move the spectrum around on the detector) were not done for these 
observations, fixed pattern noise should be present. In the other segments, some movement 
of PG 0038+199 in the LWRS aperture is expected. However, LiFlA is used for guiding, so 
the movement of the spectrum between exposures is smaller relative to the other segments. 
LiFlA data are therefore predicted to suffer the most from fixed pattern noise, hence the 
large a correction factor. 

A. 2. Correction of velocity scale 

We made trial profile fits of H2 lines with VPFIT to evaluate the accuracy of the relative 
wavelength scale within each segment. Radial velocity uncertainties for each H2 line arc 
~ 1 km s^^ (0.5 pixels). With ~ 50 H2 lines per segment (save for the long wavelength ones 
LiFlB and LiF2A), hnear dispersion corrections of 2-10 km s~^ 100 were applied, plus a 
zero point offset to make the H2 velocities match the Na I heliocentric velocity of -5.1 km s~^ 
measured from our Keck data. See Fig. 17 for an example. Standard deviations of radial 
velocity residuals after applying the linear corrections to each segment were 1-4 km s"^ 
(0.5-2 pixels). 

B. Determination of the line spread function using H2 lines 

We made an independent determination of the FUSE instrumental profile (line spread 
function, LSF), which is poorly known. The LSF is suspected to vary with detector /channel/side 
segment, wavelength, time, pointing accuracy, and how accurately individual exposures are 
cross-correlated and summed (Kruk ct al. 2002). The problem was also studied by Hebrard 
et al. (2002). Wood ct al. (2002) attempted to solve for the LSF, coming up with a double 
Gaussian fit. We attempted to incorporate their solution into the FUSE LSF for profile fits 
with VPFIT, but found that a single Gaussian gave better results based on the statistic. 
However, the FWHM of the Gaussian which gave the best fit varied between FUSE segments. 
We therefore decided to derive the FUSE LSF by using H2 lines as fiducials. This was done 
by assuming the results from a combination of profile fitting with Owens (which leaves the 
LSF as a free parameter) and the COG method, and varying the width of a Gaussian LSF 
in km s~^ for a each of a large number of H2 lines. We only used unblended lines of the H2 
J = 2,3,4 rotation levels, because the J = and 1 transitions arc in the damping part of 
the curve of growth and arc thereby mostly resolved, making them useless for assessing the 
LSF. (This is borne out by their showing uniformly infiated LSF FWHM values compared 
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to the lower column density J — 2,3,4 levels.) The number of H2 hues fitted per segment 
was 29-45 except for the two long wavelength segments (LiFlB and LiF2A), which only had 
4 and 7 respectively. 

VPFIT can accommodate a line spread function which varies as a function of wavelength, 
and uses a power series in FWHM in A as a function of wavelength. We therefore converted 
the FWHM from km to A and fitted a quadratic function to each of the segments using 
the interactive IRAF routine CURFIT, except in the case of LiFlB and LiF2A, which we 
fitted with simple constants in km s~^. The root mean square of the LSF residuals ranged 
from 0.007-0.018 A (~ 2.0 - 5.6 km s"^), with by far the worst fit for SiClB. Plots for the 
LSF determinations are in Fig. 18. 
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Table 1. PG 0038+199 stellar properties 



Quantity 



Reference 



Galactic coordinates 

Spectral type 

V 

Mass 

log(H/He) 
log(N/He) 
E(B-V) 
Distance 



119.79°, h = -42.66° 
DO 

14.544 ±0.016 
115000±11500 K 
7.5 ±0.3 
0.59 M0 
-0.7 
-3: 
<0.037 
2971}^! PC 



Green et al. (1986) 
Wesemael et al. (1993) 
Williams et al. (2001) 
Dreizler & Werner (1996) 
Dreizler Sz Werner (1996) 
Dreizler & Werner (1996) 
Dreizler et al. (1997) 
Dreizler & Werner (1996) 
Schlegel et al. (1998) 
this work 
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Table 2. H2 and HD lines used for COG analysis. 



Species 






T^lmii^;' \A/irlfli 

J— /U U.1 V V V ILi ull 


Ho 


J= 


=0 


923.990 


0.751 


125.0 ±6.3 


Ho 


J= 


=0 


931.069 


0.990 


137.7 ±13.2 


Ho 


J= 


=0 


962.985 


1.101 


200.6 ± 7.7 


Ho 


J= 


=0 


981.442 


1.306 


265.0 ±19.7 


Ho 


J= 


=0 


1077.141 


1.100 


256.8 ±8.2 


Ho 


J= 


=0 


1108.139 


0.267 


116.5 ±5.6 


Ho 


J= 


:1 


924.649 


0.580 


106.5 ±8.6 


Ho 

-^-^z 


J= 


:1 


925.180 


0.256 


98.9 ±6.9 


Ho 

^ ^z 


J= 


:1 


932.273 


0.344 


117.3 ±9.1 


Ho 


J= 


:1 


982.839 


0.824 


222.0 ±9.0 


Ho 

^^z 


J= 


:1 


992.812 


0.883 


277.8 ±11.2 


Ho 

^^z 


J= 


:1 


1003.300 


0.929 


282.6 ± 9.3 


Ho 


J= 


:1 


1108.644 


0.078 


155.9 ±6.2 


Ho 


J= 


:2 


920.248 


0.187 


38.0 ±4.2 


Ho 


J= 


:2 


927.026 


0.329 


54.9 ±4.8 


Ho 


J= 


:2 


932.614 


0.647 


61.2 ±7.8 


Ho 


J= 


:2 


933.245 


0.803 


82.6 ±6.9 


Ho 


J= 


:2 


934.151 


0.317 


45.1 ±4.1 


H2 


J= 


:2 


940.632 


0.749 


51.1 ±3.0 


H2 


J= 


:2 


941.606 


0.498 


53.7 ±3.6 


H2 


J= 


:2 


956.586 


0.941 


63.4 ±4.6 


H2 


J= 


:2 


957.660 


0.661 


58.2 ±4.6 


H2 


J= 


:2 


965.794 


1.495 


83.4 ±6.7 


H2 


J= 


:2 


967.283 


1.530 


90.9 ±5.4 


H2 


J= 


:2 


968.297 


0.844 


68.9 ±4.9 


H2 


J= 


-2 


975.351 


0.810 


62.3 ±4.1 


H2 


J= 


-2 


983.594 


1.054 


70.3 ± 3.8 


H2 


J= 


:2 


984.867 


0.905 


64.4 ± 4.6 


H2 


J= 


:2 


987.978 


1.557 


99.1 ±5.4 


H2 


J= 


:2 


989.092 


0.903 


48.2 ±4.9 



Table 2 — Continued 



Species 




loe f A^ 


Equiv Width 
fmA) 


H2 


J= 


:2 


993.552 


1.228 


96.6 ±5.9 


H2 


J= 


=2 


994.877 


0.935 


92.3 ±6.0 


H2 


J= 


=2 


1003.988 


1.222 


81.0 ±4.4 


H2 


J= 


:2 


1005.397 


0.998 


64.1 ±3.8 


H2 


J= 


-2 


1009.030 


1.200 


64.0 ±4.3 


H2 


J= 


-2 


1010.135 


1.147 


57.2 ±5.1 


H2 


J= 


-2 


1010.945 


1.381 


94.6 ±4.0 


H2 


J= 


-2 


1016.466 


1.016 


75.8 ±3.7 


H2 


J= 


-2 


1038.695 


1.234 


97.6 ± 5.2 


H2 


J= 


-2 


1040.372 


1.030 


79.5 ±4.2 


H2 


J= 


-2 


1051.505 


1.168 


77.8 ± 4.0 


H2 


J= 


-2 


1053.291 


0.980 


74.9 ± 3.7 


H2 


J= 


-2 


1065.001 


1.057 


75.8 ±4.7 


H2 


J= 


-2 


1066.907 


0.881 


85.3 ±4.3 


H2 


J= 


-2 


1079.227 


0.868 


67.9 ±5.5 


H2 


J= 


-2 


1081.268 


0.709 


79.4 ±3.7 


H2 


J= 


:3 


928.443 


0.481 


32.4 ±4.8 


Ho 


J= 


:3 


933.588 


1.260 


41.4 ±3.7 


H2 


J= 


:3 


934.800 


0.820 


51.2 ±4.8 


H2 


J= 


:3 


942.970 


0.729 


47.2 ± 3.9 


H2 


J= 


:3 


951.681 


1.079 


49.1 ±5.1 


H2 


J= 


:3 


958.953 


0.930 


50.0 ±3.9 


H2 


J= 


:3 


960.458 


0.674 


46.0 ±3.9 


H2 


J= 


:3 


966.787 


0.879 


39.5 ±5.0 


H2 


J= 


:3 


967.677 


1.340 


58.1 ±4.6 


H2 


J= 


:3 


970.565 


0.974 


53.4 ±4.5 


H2 


J= 


:3 


978.223 


0.818 


46.2 ± 3.8 


H2 


J= 


:3 


987.450 


1.406 


48.8 ±4.3 


H2 


J= 


:3 


987.772 


0.945 


51.7 ±4.9 


H2 


J= 


:3 


995.974 


1.218 


65.0 ±4.3 



Table 2 — Continued 



Species 




loe fX^ 


Equiv Width 

(uiA) 


H2 


J= 


=3 


997.830 


0.942 


59.4 ±4.9 


H2 


J= 


=3 


1006.417 


1.199 


62.4 ±4.1 


H2 


J= 


=3 


1028.991 


1.250 


55.2 ± 3.4 


H2 


J= 


=3 


1031.198 


1.059 


50.5 ±3.2 


H2 


J= 


=3 


1041.163 


1.216 


68.6 ±3.9 


H2 


J= 


=3 


1043.508 


1.052 


61.1 ±3.5 


H2 


J= 


=3 


1053.982 


1.150 


54.8 ±3.2 


H2 


J= 


=3 


1056.479 


1.006 


57.4 ±3.2 


H2 


J= 


=3 


1067.485 


1.028 


54.1 ±3.2 


H2 


J= 


=3 


1070.148 


0.909 


56.1 ±3.6 


H2 


J= 


=3 


1099.795 


0.448 


46.6 ± 3.2 


H2 


J= 


=3 


1115.907 


-0.081 


36.9 ± 3.0 


H2 


J= 


-A 


935.969 


1.264 


29.5 ±4.1 


H2 


J= 


-A 


952.765 


1.420 


33.7 ±4.1 


H2 


J= 


A 


962.158 


0.921 


19.4 ±3.6 


H2 


J= 


A 


968.671 


1.097 


26.9 ±4.3 


H2 


J= 


A 


971.392 


1.533 


33.1 ±4.1 


H2 


J= 


A 


979.808 


1.095 


22.8 ±3.4 


H2 


J= 


A 


994.234 


1.134 


21.3 ±4.4 


H2 


J= 


A 


999.272 


1.217 


23.5 ± 3.7 


H2 


J= 


A 


1011.818 


1.154 


35.8 ±3.0 


H2 


J= 


A 


1023.441 


1.031 


22.2 ±3.6 


H2 


J= 


A 


1032.355 


1.247 


21.0 ±3.0 


H2 


J= 


A 


1035.188 


1.068 


27.3 ±3.0 


H2 


J= 


A 


1044.548 


1.206 


30.1 ±3.2 


H2 


J= 


A 


1060.588 


1.019 


20.4 ±2.3 


HD J= 


-0 


959.817 


1.150 


13.4 ±4.0 


HD J= 


-0 


1007.283 


1.515 


14.8 ±3.2 


HD J= 


=0 


1011.457 


1.423 


15.9 ±2.7 


HD J= 


-0 


1042.848 


1.332 


13.6 ±2.2 



Table 2 — Continued 









Equiv Width 


Species 


A (A)^ 


log/A^ 


(mA) 


HD J=0 


1054.288 


1.238 


12.8 ±2.3 


HD J=0 


1066.271 


1.089 


6.6 ±2.1 



'^Wavelengths and /-values are from the 
Owens and vpfit atomic databases, which 
use molecular transition information from Ab- 
grall et al. (1993a,b). 
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Table 3. H2 and HD Column Densities. 



Species 


J 


Adopted (2(7)^ 


H2 





18.76 ±0.05 


H2 


1 


19.18 ±0.05 


H2 


2 


17.73 ±0.40"^ 


H2 


3 


16.74 ± 0.40^ 


H2 


4 


14.63 ±0.15 


H2 


5 


14.04 ±0.15 


H2 


total 


19.33 ± 0.04 


HD 





13.95 ±0.10 



^The values are composite results from the COG and PF 
analyses (using both Owens and vpfit). Quoted errors 
are 2(7, and we assume la errors are half the 2a uncertain- 
ties. See text for details. 



''Column density errors for H2 (J = 2, 3) are uncertain 
due to hues being on the flat part of the curve of growth. 
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Table 4. Atomic lines used for column density measurements 



S])ecies 


A (J 




log /A 


TT T 

H I 


919, 


.351 


0.043 


TT T 

H I 


920, 


.963 


0.168 


H I 


923, 


.150 


0.312 


TT T 

H I 


926, 


.226 


0.471 


TT T 

H I 


1025, 


.722 


1.909 


D I 


919, 


.102 


0.043 


D I 


920, 


.712 


0.170 


D I 


922, 


.899 


0.311 


D I 


925, 


.974 


0.470 


D I 


1025, 


.443 


1.909 


C I 


945, 


-1 r\ -1 

.191 


2.157 


C I 


1-1-1-1 

1111, 


r» -1 

.421 


0.996 


C I 


1122, 


.438 


0.820 


C I 


1129, 


-1 i~\ o 

.193 


0.998 


C I 


1139 


.793 


1.206 


C I 


1157, 


.910 


1.451 


C I 


1158, 


.324 


0.810 


C r 


945, 


.338 


2.157 


TT 

C II 


1036, 


.337 


2.102 


C III 


977, 


.020 


2.875 


TVT X 

N I 


951, 


.079 


-0.794 


N I 


952, 


.523 


-0.401 


Nil 


1083, 


.994 


2.072 


01 


974, 


.070 


-1.807 


Nal 


5891, 


.583 


3.596 


Na I 


5897, 


.558 


3.296 


Si II 


1020, 


.699 


1.188 


PII 


1152, 


.818 


2.435 


SIII 


1012, 


.501 


1.634 


Ar I 


1048, 


.220 


2.408 


Fell 


940, 


.192 


1.056 
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Table 4 — Continued 



Species 


A (A)^ 


log /A 


Fell 


1055.262 


0.926 


Fell 


1081.875 


1.180 


Fell 


1096.877 


1.545 


Fell 


1121.975 


1.351 


Fc II 


1125.448 


1.093 


Fe II 


1127.098 


0.529 


Fe II 


1133.665 


0.833 


Fe II 


1142.366 


0.757 


Fe II 


1143.226 


1.182 



'^Wavelengths and /-values 
are from the OWENS and 
VPFIT atomic databases, 
which use information from 
Morton (1991, 2003); Howk et 
al. (2000) and Verner, Barthel, 
& Tytler (1994). 
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Table 5. Atomic Column Densities. 



Ion 


Adopted 


Errors'* 


Technique^ 




Value 


J.U 






TU T 

ti I 




ZU.'4:± 


0.04 


0.08 


PF (3) 


H I+2H2 




20 48 


0.04 


0.07 


PF, COG 


F» T 
U i 




15.75 


0.04 


0.08 


Ft- y6) 


D I+HD 




15.76 


0.04 


0.08 


PF, COG 


T 




13.80 


0.04 


0.08 


Pl^ (2) 


c r 




13.25 


0.19 


+0.32 
-0.54 


PF (1) 


CII 


> 


14.47 






AOD 


cm 


> 


13.95 






AOD 


NX 




16.39 


+0.07 
-0.10 


+0.14 
-0.18 


PF (3) 


Nil 


> 


14.12 






AOD 


01 




17.37 


+0.08 
-0.15 


+0.15 
-0.25 


PF (3) 


Nal 




12.19 


0.03 


0.04 


AOD, PF (2) 


Si II 


> 


15.17 






AOD 


P II 


> 


13.34 






AOD 


S III 




13.98 


+0.03 
-0.06 


+0.07 
-0.10 


PF (1) 


Ar I 


> 


13.35 






AOD 


Fell 




14.42 


+0.03 
-0.02 


+0.06 
-0.05 


PF (3), COG 



^For comparison with other ISM work, we give 
both 1 and 2a errors. 

^PF (profile fit, Owens or vpfit, with number 
of measurements in parentheses), COG (curve of 
growth) or AOD (apparent optical depth). 
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Table 6. Unidentified lines in FUSE data for PG 0038+199 



Wavelength range (A) Wavelength range (A) Wavelength range (A) 



973.0 - 


973.2 


1062.2 - 


1062.4 


1144.7 - 


1144.9 


976.6 - 


976.8 


1067.0 - 


1067.2 


1146.0 - 


1146.5^ 


997.3 - 


997.5 


1071.9 - 


1072.1 


1146.55 


- 1146.8 


999.5 - 


999.8 


1073.65 


- 1073.85 


1147.2 - 


1147.4 


1014.6 


- 1014.8 


1080.35 


- 1080.7^ 


1147.9 - 


1148.1 


1051.9 


- 1052.1 


1080.95 


- 1081.10 


1173.45 


- 1173.89^^ 


1052.6 


-1052.8 


1125.15 


- 1125.35 


1180.48 


- 1180.68 


1061.3 


- 1061.5 


1141.1 - 


1141.3^ 







^2 hnes 

^'on the blue side of Fe II Al 144.94 
'^Fe VII Al 146.46? 
^^3 hnes 



Table 7. D I, H I, O I, N I, HD, H2 ratios. 



Species 


Ratio 


Errors (Icr) 


Errors {2a) 


10^[D I + HD]/[HI + 2H2] 


1.91 


+0.26 


- 0.24 


+0.52 


- 0.42 


lO^D/H 


2.19 


+0.30 


-0.27 


+0.55 


-0.50 


104O/[H+2H2] 


7.76 


+1.78 


-2.33 


+3.38 


-3.49 


lO^O/H 


8.51 


+1.99 


-2.56 


+4.09 


-3.86 


10^N/[H+2H2] 


0.81 


+0.17 


-0.18 


+0.36 


-0.29 


lO^N/H 


0.95 


+0.20 


-0.20 


+0.43 


-0.34 


102[DI + HD]/0I 


2.63 


+1.13 


-0.49 


+2.18 


-0.85 


lO^D/O 


2.40 


+ 1.03 


-0.45 


+2.19 


-0.78 


lO^pI + HD]/NI 


2.34 


+0.66 


-0.39 


+ 1.35 


-0.72 


lO^D/N 


2.29 


+0.65 


-0.39 


+ 1.32 


-0.71 


lO^N/O 


1.12 


+0.52 


-0.28 


+ 1.05 


-0.47 


IO6HD/H2 


4.17 






+0.68 


-0.49 
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Fig. 1.— FUSE d&tei for PG 0038+199 with ISM and geo-coronal (e) identifications. Data 
are binned by four pixels, scaled and merged for presentation only. Photospheric lines of 
Si IV are marked with dashed lines. We could not identify some absorption lines, particularly 
in the range 1080-1180 A (Table 6; a stellar origin is suspected). 
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1020 1022 1024 1026 1028 1030 
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Fig. 2. — The grid of stellar models described in § 4.1, shifted by radial velocity v = 
—52 km s^^ and plotted against the LiFlA data and la errors around H I Ly/5. Thick (dark 
grey) solid line: T = 115000 K, logg = 7.5, He-dominated (reference model). Thick (red) 
dashed lines: T = 103500 K, logg = 7.5, He-dominated (lower) and T = 126500 K, logs' = 
7.5 (upper). Thin (green) dash-dotted line: T = 115000 K, logg = 7.2, He-dominated. 
Thin (blue) dash-triple dotted line: T = 115000 K, logg = 7.8, He-dominated. Thin (light 
grey) solid line: T = 115000 K, log^ = 7.5, [H/He]=0.2, [C/H]= 5 x 10-^ [N/H]= 10"=^ 
and [0/H]= 5 x 10^^. 0: Geo-coronal LjjS, O lA 1027.5 emission. ISM and photospheric 
absorption lines are noted with ticks and labelled. The feature at 1022.7 A is not in the 
LiF2B data and appears to be noise. 
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-40 -20 20 40 

V (km s~^) 

Fig. 3. — Top panel: Apparent column density per velocity unit for the Na I doublet observed 
toward PG 0038+199. The /-values of the two transitions differ by a factor of 2. A small 
amount of unresolved saturation can be seen in the range —12 < v < —1 km s~^ where the 
profile of the stronger line falls below that of the weaker fine. Bottom two panels: Owens 
fit to the Keck HIRES data. 




Fig. 4. — Single-component COG for H2 and HD along the PG 0038+199 sight hne. A 
difTerent symbol is used for each J level. The J — 5 level was not included because the 
equivalent widths for its transitions are too small. 
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Fig. 5. — Excitation for H2 along the PG 0038+199 line of sight, using the adopted values 
from Table 3. H2 J levels are marked at top. See text for details. 



Fig. 6. — Owens profile fit to O I showing the fit before convolution with the local LSF 
(light blue line) and after convolution (red line). Flux is in erg cm~^ A"-*^. See stand-alone 
JPEG file. 
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1023 1024 1025 1026 1027 



wavelength (A) 

Fig. 7. — Solid (red) line: vpfit combined fit to H I, D I, O I and H2. Dashed (green) line: la 
error array. Dash-dotted (green) line: H2 profiles (all lines shown after convolution with local LSF). 
Dash-triple- dotted (blue) line: HI, DI profiles. Dotted (red) line: Stellar continuum. O I profiles 
are not shown because they fall in geo-coronal Ly/3 emission and are dominated by H2 lines. Q: 
Geo-coronal Ly/3, O lA 1027.5 emission. Grey zones: regions which were excluded from the fit, 
including geo-coronal emission and noise spikes. H2, HI, D I and O I transitions are marked by 
ticks. 
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Fig. 8. — Owens fit to D I Ly9 and LylO profiles (SiC2A channel, normalized data, binned 
by 3 pixels). Red lines: fit after convolution with the local LSF. Thick blue lines: individual 
components (labelled with ticks), before convolution. Solid black lines: local continua. Flux 
is in erg cm~^ s" 
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Fig. 9. — Solid (red) line: VPFIT combined fit to H I, D I and H2 (all profiles shown after 
convolution with local LSF). Dashed (green) line: la error array. Dot-dashed (green) line: 
H2 profiles. Triple- dot- dashed (blue) line: HI, Dl profiles. Dotted (red) line: Stellar contin- 
uum. Grey zones: regions which were excluded from the fit, including geo-coronal emission 
and noise spikes. The thick dashed and solid lines represent the continuum and the fit, 
respectively. H2, H I and D I transitions are marked by ticks. Flux is in units of 10~^^ erg 
cm~^ s~-^ A~^. 
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Fig. 10. — The H I 21 cm flux in units of brightness temperature vs. radial velocity toward 
PG 0038+199 from the Leiden-Dwingeloo survey. 
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Fig. 11. — Ratios, with \a and 2a errors, for [D I+HD], [H I+2H2], O I, N I column densities. 
Best fit values are shown by asterisks. 
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Fig. 12. — (a) D/H and la error bars plotted vs. line-of-sight distance, using D/H measure- 
ments in Wood et al. (2004) and this work (red). Different symbols are used for different 
sources of the D I measurement. Red open boxes: PG 0038+199, representing D/H (upper 
box) and [D I + HD]/[H I + 2H2] (lower box), slightly offset in distance for clarity, (b) D/H 
plotted vs. line-of-sight H I column density. The symbols are the same as in (a), specifically 
that the red open 6oa;es represent PG 0038-M99, showing the H I column density (upper box) 
and the [H 1+ 2H2] column density (lower box). Black horizontal and dotted lines: weighted 
means and la standard deviations as a function of A^(H I) from Wood et al. Red solid and 
dashed lines: the primordial D/H value (Kirkman et al. 2003; O'Meara et al. 2001; Pettini 
& Bowen 2001; Levshakov et al. 2002) and lower la standard deviation, respectively. For 
intermediate values 19.2 < log A^(H I) < 20.5, D/H appears variable. 
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Fig. 13. — A model curve for logA^(H I) = 20.83 plotted with the LiFlA data, which would 
be necessary to bring D/H=0.85 x 10^^, in accord with other high H I column density sight 
lines. Such a high column density is firmly excluded by the data. Line styles and symbols 
are as in Fig. 7. Column densities for D I, O I and H2 have been kept constant at their 
best-fit values. 
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Fig. 14. — DO vs. log A^(D I) for a number of sight lines, with la error bars, as in Hebrard & 
Moos (2003) and including values from Wood et al. (2004). Values for PG 0038+199 (both 
including and excluding HD contributions) are indicated by open boxes. The dotted line 
indicates the limit inside the Local Bubble, inside which the D/0 ratio is homogeneous. 



56 




13 14 15 16 

log N(DI) [cm"^] 

Fig. 15.— As in Fig. 14, but D/N vs. logA^(D I). 
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Fig. 16. — Example distribution of a error values for SiClA segment. The dotted line shows 
the point at which we set a floor for a values for SiClA data. 
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Fig. 17. — Example distribution of H2 radial velocities for segment LiF2B, which had the 
largest correction to the slope of the wavelength solution (10 km s~^ 100 A~^). The dashed 
line shows the best linear fit. The slope was removed and a linear offset applied to match 
the radial velocity derived from our Keck observations of Na I (-5.1 km s~^). 
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Fig. 18. — Fit to LSF for all segments, based on H2 lines of rotation level J = 2, 3, 4. Error 
bars are la. 
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